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It is studied to what extent the characteristics of multi-IMF (intermediate-mass 
fragments) events can be derived from the properties of the observed single-IMF 
transverse-energy spectra. It is found that the spectra of total transverse IMF en- 
ergy {Et) in multi-IMF events are well represented by "synthetic" spectra obtained by 
a multiple folding of the single-IMF transverse energy spectrum. Further, it is shown 
that, using the experimental IMF multiplicity distribution in the folding procedure, it 
is possible to reproduce the observed trends in the IMF multiplicity distributions for 
fixed values of the total transverse energy Et- Accordingly, the "synthetic" multiplic- 
ity distributions show a binomial reducibility and an Arrhenius-like scaling similar to 
that reported in the recent literature. Consistent with the concept of statistical in- 
dependence of multiple IMF production, similar results are obtained when the above 
folding-type synthesis is replaced with one based on mixing of events with different 
IMF multiplicities. For statistically independent IMF emission, the observed bino- 
mial reducibility and Arrhenius-like scaling are shown to be merely reflections of the 
shape of the single-IMF transverse-energy spectrum. Hence, a valid interpretation of 
these IMF distributions in terms of a particular production scenario has to explain 
independently the observed shape of the single-IMF Et spectrum. 

(PACS 25.70Lm, 24.10.Pa, 25.70.Mn, 25.70.Pq) 
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1. Introduction 



Recently, in several studiesBH of heavy-ion-induced nuclear multifragmentation it was con- 
cluded that certain characteristics of the intermediate-mass fragment (IMF) data are sug- 
gestive of statistical independence of multiple IMF production.! Concepts of a binomial 
or a Poissonian reducibility^ have been suggested to express the particular observed inde- 
pendence in mathematical terms. The concept of statistical independence or reducibility 
presumes that certain characteristics of multi-IMF events are reflections of the single-IMF 
production process. For example, the concept of a binomial reducibilityS presumes that the 
probability of multiple IMF production is a reflection of the single-IMF production prob- 
ability - the two quantities of interest being connected via the expression for a binomial 
distribution: 

TP ' 

n\(m — n)\ 

where P™(p) is the multiple-IMF production probability and the parameters m and p are the 
number of (binomial) tries and the probability for success in any of these trials, respectively. 

Although so far, no theoretical model has been proposed that would explain a binomial 
reducibility, one would reasonably expect that a successful model would at the same time 
link other characteristics of multi-IMF events to the characteristics of single-IMF production. 
With these expectations in mind, the present paper investigates, to what extent the empir- 
ical systematic of multiple IMF production can be derived from empirical characteristics of 
single-IMF production and, specifically, from the single-IMF transverse-energy spectrum. It 
is assumed in this study that the hypothesis of an independent IMF production holds and 
that, accordingly, the total transverse-energy spectra of multiple IMFs can be synthesized 
from the transverse-energy spectrum of a single IMF. Section 2 presents such a synthetic 
procedure, based on multiple folding of a schematic single-IMF spectrum. Section 3 dis- 
cusses results obtained when empirical single-IMF spectra are substituted for the schematic 
distributions used in Section 2, and Section 4 discusses effects of an inclusion of transverse 
energies of light charged particles in the definition of the total transverse energy. An al- 
ternative synthesis method based on fragment swapping across different events of different 
IMF multiplicity is discussed in Section 5. Section 6 summarizes the findings of the present 
study. 



2. Synthesis of Multi-Fragment E t distributions 



In the following, the total transverse energy of charged reaction products is defined as 

^ = E™f iW(©0, (2-1) 

where the summation extends over all particles with atomic numbers 1 < Z < 20. In the 
present section, only IMFs with 3 < Z < 20 are included in the definition of E t . Effects of 
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an inclusion of light charged particles are discussed further below, in Section 4. The following 
empirical characteristics of multi-IMF events represent the rationale behind the synthesis 
procedure presented in this section: 

1. The E t distributions of individual IMFs are approximately exponential, with inverse 
logarithmic-slope parameters T of the order of 30 to 50 MeVlHl As an example, the 
experimental transverse-energy spectrum is shown in Fig. 1, for single IMFs from the 
Bi+Xe reaction! at E/A = 28 MeV. Since E t is an angle- integrated quantity, the 
corresponding spectrum may differ from a Maxwell-Boltzmann shape. Superimposed 
on the data are an exponential distribution with T = 50 MeV (dotted line) and a two- 
parameter Gaussian fit (solid line) to the spectrum. While the exponential function 
deviates from the experimental spectrum both at high and low values of E t , it simplifies 
the following demonstration of qualitative features of total IMF E t spectra. 
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Fig. 1. Intermediate-mass fragment transverse-energy spectrum of single-IMF events. 
The dotted line is an exponential distribution, exp (— E t /T), used as a first-order ap- 



proximation of single-IMF spectra in Eq. The spectrum itself was used as a model 



distribution uj(e) for a numerical evaluation of the integrals ( |2.4| ) and ( |4.l| ) with a Monte 
Carlo method. 

2. To a good approximation, E t spectra of individual IMFs are independent of the IMF 
coincidence fold (multiplicity) Although discernible, the dependence on n is weak. 

3. In addition to IMFs, one also observes light charged particles (LCPs), which greatly 
outnumber IMFs, for every n. Even for n — 0, the average multiplicity of LCPs is of 
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the order of 10 and, hence, quite significant. 



4. At low bombarding energies, e.g., for E/A « 30 MeV, the IMF contribution can 
dominate the total transverse energy E t} for all but the lowest multiplicities ri.Hi'0 
For larger bombarding energies, the relative IMF contribution decreases but remains 
always significant. 

Assuming a statistically independent IMF production, the probability distribution 
Q n (ei, e n ) for events with IMF multiplicity n can be expressed in the form of a prod- 
uct of n individual single-IMF distributions Wi(ej) 

n n 

fi„(ei, ...,e n ) = n^( e ~ Il^( e *)- ( 2 - 2 ) 

t=l 1=1 

In Eq. |2.2| , as in the following discussion, single-particle variables will be denoted by 
lower-case characters, while upper-case characters indicate multi-particle variables or dis- 
tributions. For example, the symbol Q will be used to denote multi-particle probability 
distributions, while the symbol uj indicates the corresponding single-particle distribution. 
Similarly, the variable e% denotes the contribution of particle number i to the total (multi- 
particle) transverse energy E t . While the first part of Eq. |2.2| represents the statistical 
independence, the second expresses the equality of all individual energy distributions. 

Given the above form fl2.2Q of fl n (ei, e n ), the distribution fl n (E t ) = fl n (ei + ■ • ■ + e n ) 
of the total transverse energy can be expressed as 

Q n (E t ) = j "d*i ■ • • J dtn^ih, ...,t n ) ■ 5{E t - h t n ), (2.3) 

where the integration is carried out over a hyperplane of a constant E t , as enforced by 



the 5 function. Equation [O] can be further simplified to: 



r&t rti rtn-2 

ttJEt) = / oo(E t - ti)dti / u(ti - t 2 )dt 2 ■■■ v(t n -2 - t n -i)uj(t n ^i)dt n -i (2.4) 
J J o Jo 

In the general case of an arbitrary individual single-particle distribution uj, the integration 
in Eq. |2.4| can be performed numerically, at least for reasonably low IMF multiplicities. When 
a simple exponential function is taken to approximate ou, this integral can be evaluated 
analytically for arbitrarily large IMF multiplicities: 

w(ei) = ■ exp {-(3 ■ e t ) (2.5) 



n n (E t ) = (3 ■ — ■ {(3E t ) n ~ l ■ exp (-(3 ■ E t ) (2.6) 

[n — ly. 

In principle, the inverse slope parameter (3 = 1/T could depend on the IMF multiplicity 
n, but the experimental dependence on n is known to be weak. According to Eq. p^6|, 
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to the extent that the single-IMF transverse-energy spectra can be approximated by expo- 
nentials, the distribution of the composite observable E t is of Poissonian form for each IMF 
coincidence order n. 

In order to calculate quantitatively the two-dimensional joint probability distribution 
P(n, Et) for a given reaction, one should combine the one- dimensional E t distributions ( |2.6| ) 
for given n- values with the IMF multiplicity distribution y n , measured for this reaction, i.e., 

P(n,E t ) =y n -Cl n (E t ) (2.7) 

Typical experimental y n distributions are depicted in Fig. 2. For the experimental data 
discussed here, they can be approximated well by Gaussians, y n oc exp (— (n — l) 2 /2 • a 2 ) 
(solid line in the figure). 
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Fig. 2. Intermediate-mass fragment yield distribution of n-fold IMF events, for n = 0, 1, 2, 

The solid line is a Gaussian distribution, y n oc exp[— (n— l) 2 /2-er 2 ]. This yield distribution was 
used for a numerical evaluation of the integrals (|2.4|) and (JO]) with a Monte Carlo method. 



The joint distribution in Eq. |2.7| allows one to synthesize conditional IMF multiplicity 
distributions, gated with different values of E t , and to calculate the binomial parameters m 
and p for these distributions for further analysis along the lines adopted in a series of recent 
papers.0 Specifically, one calculates 



N 



N 



n 



J2i-P(hE t )/J2P(i,E t 



(2.8) 



8=1 



1=1 



N 

E(*-^) 2 
i=i 



JV 



P(i,E t )/J2 P{i,E t ), 



(2.9) 
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m = n/p, 



P = 1 - o 2 Jn. 



(2.10) 
(2.11) 
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Fig. 3. Lower panel: The distributions (2.6) are plotted up to the order of N = 11. Note 



the qualitative similarity of the distributions (2.6) and the E t distributions due to multifrag- 
mentation, plotted in the subsequent figures. Upper panel: The values of parameters p and m 
calculated with equations ( [2.11| ) and ( |2.10 ) and distributions ( |2.6| ), with sharp-cutoff weights 
y n = 1, for n up to N, and y n — for n > N. Three example values of N = 5,8, 10 are 
indicated in the figure. Note the severity of the truncation effect, which causes the calculated 
p and m to deviate from their limiting values, zero and infinity, respectively. The limits are 
reached only for N approaching infinity. 



It is the use of the constrained empirical distribution y n that leads to a reduction of the 
variances of the conditional multiplicity distributions below the Poissonian limit of n — o 2 n 
and, hence, ensures positive values for the parameters m and p. This is illustrated in Fig. 3, 
where a schematic sharp-cutoff distribution of y n serves to demonstrate the importance of 



6 



truncation. Note the large magnitude of the truncation effect, which causes the calculated 
parameters p and m to deviate from their limiting values, zero and infinity, respectively. 
These limits are reached only when the number of folds approaches infinity. Significant 
values for both p and m have been generated within the physically relevant range of E t by 
restricting the number of folds to ~ 10, i.e., a number of coincidence folds typically measured 
in an experiment. 

Expt.: weighted Poisson, T=50 MeV 
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Fig. 4. Lower panel: The E t distributions (2.7) weighted with the experimental IMF yields y n 
are compared with experimental IMF-only transverse-energy distributions. The upper part of 
the figure shows the respective p and m parameters as calculated with Eqs. ( [2. 11 ) and ( |2.1C| ), 
respectively. 



The synthetic distributions are compared in Fig. 4 with experimental IMF transverse- 
energy distributions. In the calculations, the IMF yields y n were taken from experiment .i 
The inverse slope parameter /3" 1 = 50 MeV of the IMF spectra was adjusted to approx- 
imately match the singles (n = 1) spectrum up to about 200 MeV (see Fig. 1). As seen 
in this figure, already a synthesis using an exponential singles spectrum reproduces quali- 
tatively the gross features of the experimental E t distributions. The lack of a quantitative 
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agreement with experimental data is not surprising, since the used schematic exponential 
single-particle distribution ( |2.5| ) is only a very crude approximation of the actual single- 
particle spectrum. As seen in Fig. 1, the tail of the exponential function does not match the 
n = 1 spectrum for E t > 200 MeV. Nevertheles s, the overall features of the experimental 
E t distributions are reproduced qualitatively by those calculated with Eq. pTT[ 

The upper part of Fig. 4 presents a comparison between the parameters p and m, ex- 
tracted from experimental data, with the ones calculated from Eqs. [2.8| and |2.9| , using the 
distributions (|2.7| ) displayed in Fig. 2. The predicted values of p and m match the experi- 
ment remarkably well, even though the calculated E t distributions in the lower panel deviate 
from the experimental ones. This points towards the need for a better description of exper- 
imental single-fragment E t distributions u(e). More realistic shapes of these single-particle 
distributions will be considered next. 

3. Synthesis with Realistic Single IMF Distributions 

A more realistic rendering of the shape of the single-IMF distribution than achieved by an 
exponential is provided by a two-parameter Gaussian fit of the kind depicted by the solid line 
in Fig. 1. Such a Gaussian spectrum was used in a numerical evaluation of the integrals Q2.4| ), 
employing a Monte-Carlo sampling method results of which are displayed in Fig. 5. As seen 
in Fig. 5, a remarkable quantitative agreement with experimental data is achieved. This 
time, also the parameters p and m have been reproduced. It is worth emphasizing, that 
only four parameters, two shape parameters of the single-particle spectrum (cf. Fig. 1) and 
two parameters describing the IMF yield distributions (cf. Fig. 2)) are sufficient to reproduce 
the shapes of the multi-IMF E t distributions, and the E^-dependencies of the parameters p 
and m. 

4. The Role of LCP Contributions 

In order to better represent the experimental situation, the contribution by light charged 
particles to the total kinetic energy has to be considered. This can be achieved easily 
by extending the concept of a statistically independent fragment production to the LCP 
emission. In this case, 

a n {E t ) = tt{e + e x +, +e n ) = uj {E t - t )dt f° co(t - t^dtx.... (4.1) 

Jo Jo 

The index now denotes the LCP contribution to the total E t . Since there are many 
LCPs contributing independently (even for the case of no IMFs, n — 0), according to the 
Central Limit Theorem, the contribution of all LCP's, taken together, is approximately a 
Gaussian 

uj (e ) oc exp [-(e - e^) 2 /2a 2 }, (4.2) 



S 



where the parameters eo and ctq may depend on the IMF multiplicity n. 



Expt.: single IMF Gaussian fold 
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Fig. 5. Lower panel: IMF transverse- energy distributions (solid lines), predicted by folding 
a Gaussian form of the n — 1 IMF distribution multiple times with itself (without the LCP 



contribution), according to formula (2.4). The resulting E t spectra were weighted with the 
Gaussian IMF yields y n from Fig. 2. The symbols connected with dotted lines represent 
experimental data. Upper panel: p and m parameters extracted from the experimental data 
(symbols connected with lines), compared with the ones extracted from the set of folded E t 
distributions. 



Wieloch et. al.u'Q have considered the particular parameterization of e5" oc n and do = 
const. The present paper employs values of eo and a Q taken from experiment ,0 namely 
eo oc n for n = 0,1,2 and eo ~ const for n > 2. With this choice of parameter values 
eo and ctq, the integral ( |4.1|) was evaluated numerically in the same way as before. The 
results are displayed in Fig. 6, showing satisfactory agreement with experimental data. 
The parameters p and m are also rather well reproduced by the calculation, except for 
the deviations from experimental data at E t < 100 MeV. The deviations are attributed to 
an imperfect representation of LCP spectra by Gaussians, and in particular to the intense 
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contributions of peripheral events to the LCP spectra at low E t . It should be noted, that 
low-Et data are usually considered less important in the context of multifragmentation. 



Synth.: n IMF 
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Fig. 6. Lower panel: IMF transverse- energy distributions (solid lines), predicted by folding a 
Gaussian approximation to the E t distribution for n = 1 multiple times with itself and with the 



LCP contributions, according to formula (4.1). The resulting E t spectra were weighted with 
the yields y n from Fig. 3. The symbols connected with dotted lines represent experimental 
data. Upper panel: Experimental p and m parameters are compared with those predicted by 
folding. 



5. Mixed-event calculations 

As an alternative approach to synthesizing multi-IMF events, the experimental data! have 
been subjected to "event mixing", a procedure where both IMFs and LCPs are randomly 
exchanged between different events. For every event, the number of LCPs, ulcp, and the 
number of IMFs, n, were preserved, but not the identities of these particles. Both IMFs and 
LCPs were swapped among events of all multiplicities (regarded as a common "fragment 
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pool"), in order to ensure a complete "loss of memory" in the synthetic mixed-event data 
set. 

The above random particle exchange between events ensured that multi-particle corre- 
lations were erased from experimental data, while preserving approximately the shape of 
the single-particle distributions. In this fashion, a statistical independence among particles 
was strictly imposed on experimental data, even if some degree of multi-particle correlations 
had been present in the original data set. This procedure ensured that the assumptions 



underlying Eq. |4.1| were strictly fulfilled by the data. It should be stressed, that event mix- 
ing does not represent any additional assumption brought into the problem. Rather, it is 
a convenient way to evaluate the integral (|4.1| ), where all single-particle probability distri- 
butions were sampled directly from the experimental data, without recourse to a particular 
parameterization. 

Synthetic: mixed events 
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Fig. 7. Transverse-energy distributions (lower panel) and p and m parameters (symbols with 
dotted line in upper panel) extracted from the experimental data set and from the same set, 
but subjected to random particle swapping (solid lines). 



Subsequently, the synthetic events were subjected to the same analysis procedure as the 
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original data. The results of this mixed-event analysis are compared with the original data 
in Figs. 7 and 8. Very good agreement between both sets of results is observed from these 
figures, both with regard to the E t distributions and the extracted parameters p and m. 
This result confirms again that the concept of a statistically independent IMF production 
is a valid one. 
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Fig. 8. The p and m parameters extracted from the experimental data set (symbols connected 
with solid line) and from the same set subjected to random particle swapping (dotted lines), 
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plotted in an Arrhenius-like representation as functions of E t . 




6. Summary 

Accepting the validity of the findings of earlier papersB! regarding the statistical indepen- 
dence, or reducibility, of multiple intermediate-mass fragment production, one finds that 
important characteristics of multi-IMF events can be derived already from the experimental 
shapes of single-IMF transverse energy distributions. The characteristics of the synthetic 
events generated from these distributions show remarkable resemblance to those of the ex- 
perimental multi-IMF events. The latter characteristics include the binomial reducibility 
and an Arrhenius-like scaling discovered in recent studies of nuclear multifragmentation.0 

Although the findings of the present study do not directly favor one IMF production 
scenario over another, they do so indirectly. They draw attention to the crucial importance 
of an understanding of the shapes of the single-IMF E t distributions. In the present analysis, 
the apparent binomial reducibility and the apparent Arrhenius-like scaling are secondary 
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phenomena, both reflecting the "input" shape of the elementary E t distribution. Hence, a 
valid interpretation of the significance of such reducibility and scaling in terms of a particular 
IMF production scenario requires an independent explanation of the underlying single-IMF 
E t distribution. 

A purely thermal scenario, which could potentially justify Arrhenius-like correlations 
between emission probability and temperature of the system, makes quite definite predic- 
tions regarding the elementary E t distributions, which should be very similar to those of 
the lighter particles emitted from the same sources. These predictions disagree with the 
experimental IMF spectra, which are significantly harder than the light-particle E t spectra. 
This observation justifies some skepticism towards an interpretation!! of Arrhenius-like plots 
in terms of thermal scaling. 
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